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Abstract 

In the papers IS El we have established linear and quadratic bounds, in k, on 
the growth of the Lebesgue constants associated with the /c-sections of Leja sequences 
on the unit disc U and JI-Leja sequences obtained from the latter by projection into 
[—1,1], In this paper, we improve these bounds and derive sub-linear and sub-quadratic 
bounds. The main novelty is the introduction of a “quadratic” Lebesgue function for 
Leja sequences on U which exploits perfectly the binary structure of such sequences and 
can be sharply bounded. This yields new bounds on the Lebesgue constants of such 
sequences, that are almost of order y/k when k has a sparse binary expansion. It also 
yields an improvement on the Lebesgue constants associated with J?-Leja sequences. 


1 Introduction 

The growth of the Lebesgue constant of Leja sequences on the unit disc and jft-Leja sequences 
was first studied in 011]. The main motivation was the study of the stability of Lagrange 
interpolation in multi-dimension based on intertwining of block unisolvent arrays. Such 
sequences, more particularly jft-Leja sequences, were also considered in many other works 
in the framework of structured hierarchical interpolation in high dimension. Although not 
always referred to as such, they are typically considered in the framework of sparse grids 
for interpolation and quadrature mm- Indeed, the sections of length 2 n + 1 of jft-Leja 
sequences coincide with the Clenshaw-Curtis abscissas cos(2 -n j7r), j = 0,..., 2 n which are 
de facto used, thanks to the logarithmic growth of their Lebesgue constant. 
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Motivated by the development of cheap and stable non-intrusive methods for the treat¬ 
ment of parametric PDEs in high dimension, we have also used these sequences in 5j 
with a highly sparse hierarchical polynomial interpolation procedure. The multi-variate in¬ 
terpolation process is based on the Smolyak formula and the sampling set is progressively 
enriched in a structured infinite grid <8>j =0 Z together with the polynomial space by only one 
element at a time. The Lebesgue constant that quantifies the stability of the interpolation 
process depends naturally on the sequence Z. We have shown in [?] that it has quadratic 
and cubic bounds in the number of points of interpolation when Z is a Leja sequence on 
U. or an SK-Leja sequence, thanks to the linear and quadratic bounds on the growth of the 
Lebesugue constant of such sequences, also established in [6j[?]. We refer to the introduction 
and section 2 in [7j for a concise description on the construction of the interpolation process 
and the study of its stability. 

The present paper is also concerned with the growth of the Lebesgue constant of Leja 
and SK-Leja sequences. We improve the linear and quadratic bounds obtained in [7J. In 
particular, we show that for JNLeja sequnences, the bound is logarithmic for many values of 
k which may be useful for proposing cheap and stable interpolation scheme in the framework 
of sparse grids tm 

1.1 One dimensional hierarchical interpolation 

Let X be a compact domain in C or R, typically the complex unit disc U := {|z| < l}or the 
unit interval [—1,1], and Z = o a sequence of mutually distinct points in A". We denote 
by Iz k the univariate interpolation operator onto the polynomials space Pfc_i associated with 
the section of length k, Z & := (z 0 , • • •, Zk- 1 ). The interpolation operator is given by Lagrange 
interpolation formula: for / G C(X) and z € X 

k—1 k —1 

I zJ(z) = '52f{z j )l jik {z), lj,k(z ) =Y[ Z Z \ j = 0 , • • •, k - 1 . ( 1 . 1 ) 

j =0 i=0 ** z 3 

Since the sections Z^ are nested, it is convenient to give the operator Iz k using Newton 
interpolation formula which amounts essentially writing: A 0 / := IzJ = f{z o) and 

k 

Iz k+1 = Iz k + Afc = ^ A I, where A t (Z) := I Zl+1 — Iz v l > 1, (l- 2 ) 

1=0 

and computing the operators A i using divided differences, see [Oj Chapter 2] or equivalently 
the following formula which are differently normalized, see Hi, 

a,/=(/ w-fe I /^))n[ £ A, i>i. (i.3) 

v J JIq l Z 1 - z i) 
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The stability of the operator I Zk depends on the positions of the elements of Z k on A", 
in particular through the Lebesgue constant associated with Z k , defined by 


L z k ■ = max 

feC(x)-{ 0 } 


\\IzJ\\ 


L°°(X) 


L°°(X) 


= max A Zk (z), 


zex 


where X Zk is the so-called Lebesgue function associated with Z k defined by 


(1.4) 


fc-i 

X z k {z) := ^2\h,k(z)l zeX. 

i =0 


We also introduce the notation 

B k (Z) 


max 

feC(x)-{ o} 


Afc/|U°°(x) 

ll/IU°°(X) 


(1.5) 


( 1 . 6 ) 


In the case of the unit disk or the unit interval, it is known that L*. the Lebesgue constant 
associated with any set of k mutually distinct points can not grow slower than 21 °]^ and it is 
well known that such growth is fulfilled by the set of fc-roots of unity in the case X — U and 
the Tchybeshev or Gauss-Lobatto abscissas in the case X = [—1,1], see e.g. [2], However 
such sets of points are not the sections of a fixed sequence Z. 

In [3 3], the authors considered for Z Leja sequences on U initiated at the boundary 
dU and 3?-Leja sequences obtained by projection onto [1,1] of the latter when initiated 
at 1. They showed that the growth of h Zk is controlled by 0{k,\og{k)) and 0(k 3 log(fc)) 
respectively. In our previous works M, we have improved these bounds to 2 k and 8y/2k 2 
respectively. We have also established in |7] the bound O*, < (1 + k) 2 for the difference 
operator, which could not be obtained directly from O fc < L Zk+1 + h Zk and which is essential 
to prove that the multivariate interpolation operator using 3ft-Leja sequences has a cubic 
Lebesgue constant, see [7, formula 25]. 


1.2 Contributions of the paper 

In this paper, we improve the bounds of the previous paper I3E1IS1E ]. Our techniques 
of proof share several points with those developed in [M, yet they are shorter and relics 
notably on the binary pattern of Leja sequences on the unit disk. The novelty in the present 
paper is the introduction of the “quadractic” Lebesgue constant 

k —1 i 

^z k ,2(z) ■= \h,k{z)\ 2 ^j , zex, (1.7) 

i =0 

where l l>k are the Lagrange polynomials as defined in fll.51) . We study this function and its 
maximum 

h Zk ,2 := maxAz fe ,2(^)- (1-8) 

z£X 
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We establish in §2 in the case where Z is any Leja sequence on U initiated on the boundary 
dU the “sharp” inequality 

A z k {zk ) < ^z k ,2 < 3Az k (z/c), arid A z k {zk) '■= V2 a X k ) — 1 ; (1.9) 

where (Ji(k) denote the number of ones in the binary expansion of k. Cauchy-schwatrz 
inequality applied to the Lebesgue function X Zk defined in (jl.5j) yields A z k < Vk Xz k , 2 - This 
shows that we also establish 

h Zk < 3 Vk v 7 2 a X k ) - 1, (1.10) 

for Leja sequences on U, which improves considerably the linear bound 2k established in [6] 
when the binary expansion of k is very sparse. For example, for k = 2 n + 3 with n large, 
we get h Zk < 3\/7 k « k. Using the bound (jl.9j) . we establish in §3 a new bound on the 
growth of Lebesgue constants of K-Leja sequences that implies 

L fc < 6^5 k 2 <TlW ; w h er e l = k - (2 n + 1), (1.11) 

where n is the integer such that 2 n + 1 < k < 2 n+1 + 1. Again, we remark that the previous 
bound improves the bound 8x/2. k 2 established in [7j when l is small compared to 2 n or very 
sparse in the sense of binary expansion. We actually prove a bound that is logarithmic for 
many values of k other than the values 2 n + 1, see Theorem 13.21 

Finally, we provide in 0 new bounds on the growth of O*. the norm of the difference 
operators. We provide the bounds 

D fc < 1 + ^k(2 a Xk) _ i) ; D fe < 2 <Tl(fc) 2 n , k> 1. (1.12) 

in the case of Leja sequences on U and the case of K-Leja sequences respectively where for 
the latter n is defined as above. 

1.3 Notation 

In the remainder of the paper, we work with the following notation. For an infinite sequence 
Z \= (zj)j >o on X, we introduce the section Zi m := (zi, ■ ■ ■, z m _i) for any l < m — 1. 
Given two finite sequences A = (a 0 ,..., a^- 1 ) and B = (6 0 , • • •, h-i), we denote by A A B 
the concatenation of A and B , i.e. A A B = (do,..., ak-i, b 0 ,..., 6;_i). For any finite set 
S = (so, • • •, Sfc-i) of complex numbers and p G C, we introduce the notation 

pS := (ps 0 , • • • ,psfc_i), K(S') := (3?(s 0 ),---,3fJ(s fc _i)), S := (sjf, • • •, sTm). (1-13) 

Throughout this paper, to any finite set S of numbers, we associate the polynomial 

Ws(z ) := J^[(^ — s), with the convention w$(z) := 1 (1.14) 

s£S 
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Any integer k > 1 can be uniquely expanded according to 

n 

k = ^a j 2 j , dj G {0,1} (1.15) 

3=0 

We denote by (k), ao (k) the number of ones and zeros in the binary expansion of k and 
by p(k ) the largest integer p such that 2 P divide k. For k = 2 n ,..., 2 n+1 — 1 with binary 
expansion as above, one has 

n n 

a i (k) = ^^aj and a 0 (k) = ^^(1 — dj) = n + 1 — ^(k). (1.16) 

3=0 j=o 

We should finally note that, unless stated otherwise, we only work with complex numbers 
z belonging to the unit circle dU. This is because in the complex setting we investigate supre- 
mums of sub-harmonic functions, A z k and \z k , 2 i which is always attained on the boundary. 

2 Leja sequences on the unit disk 

Leja sequences E = {ej)j>o on U considered in [3| 0, [ 6 ], 7] have all their initial value eo G dU 
the unit circle. They are defined inductively by picking eo G dU arbitrary and defining e k 
for k > 1 by 

e k = argmax^Jz - e fc _i|... \z - e 0 \. ( 2 . 1 ) 

The maximum principle implies that ej G dU for any j > 1. Also, the previous argrnax 
problem might admit many solutions and e k is one of them. We call a fc-Leja section every 
finite sequence (e 0 ,..., e k -i) obtained by the same recursive procedure. In particular, when 
E := (ej)j>i is a Leja sequence then the section E k = (eo,..., e k -i) is fc-Leja section. 

In contrast to the interval [—1,1] where Leja sequences cannot be computed explicitly, 
Leja sequences on dU are much easier to compute. For instance, if eo = 1 then we can 
immediately check that e\ = — 1 and Assuming that e 2 = i then e 3 maximizes 

\z 2 — l\\z — i\, so that 63 = —i because —i maximizes jointly |z 2 — 1| and |z — i\. Then 
e 4 maximizes |z 4 — 1|, etc. We observe a “binary patten” on the distribution of the first 
elements of E. 

By radial invariance, an arbitrary Leja sequence E = (e 0 ,ei,...) on U with e 0 G dU is 
merely the product by eo of a Leja sequence with initial value 1. The latter are completely 
determined according to the following theorem, see PI 

Theorem 2.1 Let n > 0, 2 n < k < 2 n+1 and l = k — 2 n . The sequence E k = (eo,..., e^-i), 
with eo = 1, is a k-Leja section if and only if E 2 n = (eo,..., e 2 «-i) and Ui = (e 2 »,..., efc-i) 
are respectively 2 n -Leja and l-Leja sections and e 2 ™ is any 2 n -root 0 /— 1. 
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In the light of the previous theorem, a natural construction of a Leja sequence E := (ej)j >o 
in U follows by the recursion 

E\ := (eo = 1) and E 2 n+i := E 2 n A e^^", n > 0. (2-2) 

This recursive construction of the sequence E yields an interesting distribution of its ele¬ 
ments. Indeed, by an immediate induction, see [I], it can be shown that the elements e k are 
given by 

n n 

e k = exp ^Z7r a j2 ~ :i j for k = ^ ap? ■, aj G {0,1}. (2.3) 

j=o j=o 

The construction yields then a low-discrepancy sequence on dU based on the bit-reversal 
Van der Corput enumeration. 

As already mentioned above, Theorem 12.11 characterizes completely Leja sequences on 
the unit circle. It has also many implications that turn out to be very useful in the analysis 
of the growth of Lebesgue constants. 

Theorem 2.2 Let E : = {ej)j> 0 be a Leja sequence on U initiated at e 0 G dU. We have: 

• For any n > 0, E 2 n = e 0 U 2 n in the set sense where U 2 n is the set of 2” -root of unity. 

• For any k> 1, \w Ek (e k )\ = sup zedu \w Ek (z)\ = 2 ffl(fc) . 

• For any n > 0, E 2 n, 2 n+i := (e 2 ™, • • •, e 2 ™+i_i) is a 2 n -Leja section. 

• The sequence E 2 := (e 2 j)j>o is a Leja sequence on dU. 

Such properties can be easily checked for the simple sequence defined in (12.311 and are 
given in M for more general Leja sequences. 

2.1 Analysis of the quadratic Lebesgue function 

It is proved in [6] that given two k- Leja sections Ej, and Fk, one has F-,\ = pEk in the set 
sense for some p G dU. This means that the sequence F k can be obtained from E k by a 
permutation and the product by p. By inspection of the quadratic Lebesgue function (II.8p . 
we have then that 


^F k ,2(z) - \ Ek ,2(z/p), z eU => h Fk)2 - h Ekt2 . (2.4) 

In order to compute the growth of Lg fci2 for arbitrary Leja sequences E, it suffices then to 
consider E to be the simple sequence given by (12.31) . Unless stated otherwise, for the rest of 
this section, E is exclusively used for this notation. Let us note that 

E 2 := (e 2 j)j> o = E. (2.5) 
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In order to study the functions \E k , 2 , we adopt the methodology that we introduced in 
| 6 j. Namely, we study the implication of E being a Leja sequence in general, on the growth 
of A 7 ^ 2 , then we use the implication of the particular binary distribution of E to derive such 
growth. 

Lemma 2.3 Let Z be a Leja sequence on a real or complex compact X. For any k > 1 and 
any z G X, it holds 

^z k+1 p(z) < A Z k ,2( z ) + A Z k g{ z k ) + 1- (2-6) 


Proof: We fix k > 1 and denote by Iq, ..., l k ~i the Lagrange polynomials associated with 
the section Z k and by L 0 ,... ,L k the Lagrange polynomials associated with the section Z k+l . 
By Lagrange interpolation formula, for j — 0,..., k — 1 


k 

lj ^ ^ l ] {z i )L l Lj + lj ( z k j L k 
i =o 


We have then for any z G X 


lj lj i z k)L k . 


Qr i Ljwr) 172 < (f) i/,(*)i 2 ) 1/2 +ii*wi(£ i^-mi 


k -1 


fc -1 


1/2 


j =0 


3=0 


j =0 


where we have merely applied triangular inequality with the euclidean norm in C fc . This 
also writes 


|Az fc+1 , 2 ( 2)1 - \L k (z)\ J < A Zfcj2 ( 2 ) + \L k (z)\\ Zk p(z k ). 

We conclude the proof using a < \Ja 2 — b 2 + b for a > b > 0, and the inequality 

I t (A \ _ \ w z k {z)\ 

\l J k\ z )\ — I / 11 S f 

\wz k (z k ) I 


which follows from the Leja dehnition (12.11) . 


The previous result shows that given Z a Leja sequence over X , the growth of L ^ fc)2 is 
monitored by the growth of A z k , 2 (z k ). In particular, it is easily checked using induction on 


k that 

Az fc , 2 ( 2 *;) = 0(log(k)) = 

=► L z , j2 = 0(kbg(k)), 

(2.7) 

and 

A zM = 0{k e ) = 

=> u h ,i = 0(V+'). 

(2.8) 
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In the following, we show basically that the previous implication holds with 6 = 1/2 for 
Leja sequences on U. However, we use the particular structure of such sequences in order to 
show that the exponent 6 = 1/2 is not deteriorated and that it is also valid for Lg fc) 2 . We 
recall that we work with the simple sequence E given in (12.31) for which E 2 = E. The binary 
patten of the distribution of E on the unit disc yields the following result. 

Lemma 2.4 Let E be as in (12.31) . For any N > 1, one has 

^e 2N ,2(z) = ^e n ,2{z 2 ), z G dlA. (2-9) 


Proof: Let Iq, ..., I 2 N -1 be the Lagrange polynomials associated with E 2 n and L 0l ..., Lat_ 1 
be the Lagrange polynomials associated with E N . Since e 2j+ i — —e 2 j for any j > 0, then in 
view of (j2.5[) 

we 2 N (z) = w E 2 f (z 2 ) = w En (z 2 ). 

Deriving with respect to z and using (e 2 j + i) 2 = (e 2 j) 2 = for any j > 0, we deduce that 


Wt 


J, E 2N ( e 2j+i)l = \ w E 2N ( e 2j)\ = 2| w' EN {e\j)\ = 2|^(ej)|, j > 0. 
We have for any j = 0,..., N — 1 

\we 2N {z)\ 1 , _ \we 2N (z)\ 


( 2 . 10 ) 


1^2 j (^) | | / / \ 11 | j 

\ W E 2N ( e 2j)\\ Z - e 2jl 

Therefore in view of the previous equalities 


\hj+i(z)\ — 


WE 2N ( e 2 j+ l)\\z - e 2j+l\ 


2\\2 r 


\l 2j {z)\ 2 +\l 23+l {z)\ 2 = 


M z2 )\ 2 


\we n (z 2 )\ 

M w E N ( e j) I 2 L k - e 2i| 2 ' \z + e 2j \ 2 \ We {ej)\ 2 \z 2 - e 


= I LAz 


2\ 12 


( 2 . 11 ) 

where we have used \a — b\ 2 + \a + b\ 2 = 4 for a, b G dU and e 2 ^ = e r Summing the previous 
identities for the indices j — 0,..., N — 1, we get the result. ■ 


We note that the previous result combined with E 2 n = U 2 n in the set sense implies that 


2 n — l 


E 


2 n 1 

z J —1 

2 n (z - 6j) 


A E 2 n,2(z ) — \e!,2{z 2 ) — 1, 


( 2 . 12 ) 


for any z G dlA. We now turn to the growth of \E k , 2 (ek), which as mentioned earlier monitor 
the growth of L Efei2 . 
















Lemma 2.5 For the Leja sequence E defined in (12.31) . we have for any k > 1 


A E .,2(eib) = -1 


(2.13) 


Proof: First, by Lemma [2.41 and e\ N — ejv, one has 

|A_B 2JV ,2(e2Ar)| 2 = |AEjv,2(eAr)| 2 , iV > 1. 


(2.14) 


Let now k be an odd number and we write k = 2N + 1 with N > 1. Let l 0 ,..., 1 2 n be the 
Lagrange polynomials associated with E k and L 0 ,..., L N -\ be the Lagrange polynomials 
associated with E^. For any m = 0,..., 2N, one has 




w Ek {ek) 


w 


E k +i 


(e k ) 


(e k - e m) w E k (e m ) w' Ek+i (e m ) 


I Im ( e k ) | 


Wt 


^N+l 


[el) I 


\w 


En+i 


[<)V 


where we have used k + 1 = 2(N + 1) and (12.101) . Using e\ = e n and (e 2 j+i) 2 = (e 2 j) 2 = Cj 
for any j, we get for m = 2j or m = 2j + 1 with j — 0,..., N — 1 


|^m(Cfc)| 


\W 


En+i 


( e iv)l 


to 


^N +1 


(e S )\ 


= | Lj(e N )\ and also |Z 2 Jv(efc)| = 


Wt 


,( e Ar)| 


\Wi 


^N +1 


(ew)| 


= 1. 


Summing the numbers \l m (e k )\ 2 over m = 0,..., 2N , we infer 

|AE 2 jv+i,2( e 2V+i)| 2 = 2|AE iV ,2(ev)| 2 + 1- (2.15) 

In view of the above and A^ li 2 (ei)=l, the sequence a := ( a k := |A^^(e*,)| 2 )fc>i satisfies: 
ex i = 1 and ct 2 N = c^N: & 2 N +1 = 2ctjv 4*1; IV > 1. 


We have Oi(l) = 1 and cr 1 (2iV) = a\(N), ai(2N + 1) = cri(iV) + 1 for any N > 1. It is 
then easily checked that (2 CT1 ^ fc ' ) — l)fc>i satisfies the same recursion as a. This shows that 
a k = 2 CT d fc ) _ l for any k > 1 and finishes the proof. ■ 


We are now able to conclude the main result of this section, which states basically that 
for the sequence E or more generally any Leja sequence on U initiated at the boundary dU , 
the value of L E k ,2 = max ze ^ \E k , 2 (z) is almost equal to A E k , 2 (e k )- 


Theorem 2.6 For the Leja sequence E defined in (12.31) . we have for any k > 1 

^ ^ ^E k ,2 _ jUgfc, 2 ^ 0 

^E k , 2(e k ) y/2 a F k ) — 1 


(2.16) 
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Proof: The first part of the inequality is immediate from the definition of L E k , 2 - Also in 
view Lemma 12.41 and formula (12.14)1 . we only need to show (12.1611 when k is an odd number. 
Let k = 2N + 1 with N > 1. Using Lemma [2.31 Lemma [2.41 and formula (12.14[) . we have 

^E k , 2 ( 2 ) < ^E 2N , 2 (A) + ^ E 2N ,‘2( e ‘2N ) + 1 = A^ JV]2 (^ 2 ) + A E N ,2( e N ) + 1- 
If we assumes that Xe n ; 2 (z 2 ) < 3A e n , 2 (zn), we get 

2 ( 2 ) < dAe^^e^) + 1 < 3-^2|A^;^^(e^v)| 2 + 1, 

where we have used the elementary inequality At + 1 < 3\/2 1 2 + 1 for any t > 0. In view 
of (I2.15jl . one then gets A E k V z ) < 3XE k ,2( e V- The verihcation L# 1)2 = Ae li2 (ei) = 1 shows 
that the result follows using an induction on k > 1. ■ 


2.2 Implications on the Lebesgue constant 

The methodology we have provided so far for bounding L E k) 2 is not new, we have developed 
it in |6j in order to give linear estimate for L E k , namely L E k < 2 k. Theorem (I2.6jl has 
also implications on the growth of the Lebesgue constant L E k - Indeed, Cauchy Schwartz 
inequality applied to the Lebesgue function X Ek implies A< Vk A E k , 2 i so that 

h Ek < Vk h Ek , 2 < 3^(2-d fc ) - 1). (2.17) 

The Cauchy Schwartz formula \ Ek < Vk A E k2 is possibly not very pessimistic. It has been 
recently proved that the Lagrange polynomials are uniformly bounded, see [Qj We shall 
observe in particular, see Figure, that the binary pattern observed for the exact value of L E k 
is captured by the previous bound. Moreover, we are able to provide a lower bound for L E k , 
that is comparable to the previous upper bound for values of k with full binary expansion. 


Proposition 2.7 For the Leja sequence E defined in (j2.3f) . we have for any k > 1 


- 1 < \ Ek (e k ) < h Et . 


(2.18) 


Proof: We let IV > 1 and we use the notation of the proof of Lemma 12.41 As for formula 
(12.lip and since |a — b\ + \a + b\ >2 for any a, b e dU , one has 


\hj(z)\ + \hj+i(z)\ 


\we n (z 2 )\ \z-e 2j \ + \z + e 2j \ > , 2 

2 l w W( e i)l \ z ~ e j\ 


This implies A e 2 N (z) > A e n (z 2 ) and more particularly A e 2N (s 2 n) > A e n Vn)- As in the proof 
of Lemma [231 we have also A£ 2JV+1 (e 2JV+ i) = 2A Ejv (e iV ) + L The sequence (b k \E k (e k )) k >i 
satisfies: 


b\ — 1 and & 2 v ^ 6 jv ? ^ 2 ./v+i — 2+ 1, AC > 1. 
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The sequence b then satisfies > 2 CT d fc ) — 1 for any k > 1. 


The previous theorem combined with Theorem 12.61 and (12.17ft implies 

1 /2 cri (*0 _ l 

-2- ^E kl 2 < I -E k < Vk L Efej2 . (2.19) 

Cauchy Schwartz inequality is then satisfactory when k ~ 2 ai ^ k \ that is when k has a full 
binary expansion. 

Remark 2.8 For integers k = 2 n ,..., 2 n+1 — 1, if k = 2 n+l — 1 in which case a\(k) — n + 1 
is the largest possible, the bound (12.17li merely implies L E k < 3 k which is worse than the 
bound 2k established If and the exact value h Ek = k of this case, see f3j/. However, since 
ai(fc) = n + 1 — ao(k) for any k = 2 n ,... 2 n+1 — 1, then by (I2.17P 

I 18 I 18 

< 2 - 20 > 

This shows in particular that L Ek < k whenever cr 0 (k) > 5. This last result answers partly 
the conjecture raised in m and which states that h Ek < k for any k > 1. 

For the purpose of the next section, we improve the bound (j2. 1 7j) in the case where k is 
an even number. We recall that we have shown in |[6J Theorem 2.8] 

^e 2 pi < h 2P h El , p > 0, l > 1, (2-21) 

where L 2P is the Lebesgue constant associated with the set of 2 p -roots of unity. The value 
L 2p can be computed easily for small values of p and it grows logarithmically in 2 P , see e.g. 
|6j formula 2.25], 

L x = 1, L 2 = y/2, and L 2P < l(log(2 p ) + 9/4), p > 2. (2.22) 

Since a\(k) = ai{k/2 p ^), we have then in view of (12.171) and (I2.2ip the following theorem 

Theorem 2.9 Let E be the Leja sequence defined in (12.3ft or any Leja sequence onU initiated 
at dU. We have 

Le„ < 3yC|_(2».(‘> - 1) L 2 „ w , k> 1. (2.23) 


We should mention that our primary interest in studying X Ek} 2 was the improvement of 
the results of [7] concerned with the Lebesgue constants of K-Leja sequences. This will be 


11 



























made clear in the proof of Theroem 13.21 For the sake of the same theorem, we need also to 
provide a growth property of Leja sequences on the unit disc. 

We let E = (ej)j>o be the simple Leja sequence defined by (12.3[) . For m > 0 and 
1 < l < 2 m_1 , we introduce the notation K = 2 m + l and F m j = _F 2 m ,A' and define the 


quantity 


'Ym,l . m 


K -1 

E 

3 =0 



(2.24) 


The quantity 7 m ,z is well defined. Indeed, by the particular structure of the sequence E : we 
have E- 2 m, + 2 m-i = E 2 m so that Fi 2 m + 2 m - 1 = W 2 ™ Ae^ U 2 m -1 in the set sense. We 

have then for j = 0,..., 2 m + / — 1, ej is in W 2 m A e 2 ^ U 2 m - 1 which does not intersect with 
F m j C U 2 m-i. We have the following growth for y m j. 


Lemma 2.10 For any m > 1 and any 1 < l < 2 m we /love 

. 5 

7m ,1 — 2 < t 1 ( 0 + p ( 0 +i 


(2.25) 


Proof: Since (e 0 ,ei,e 2 ) = (1, — 1, z), it can be checked that 7^1 = 5/4. We then fix m > 2. 
We define p = e 2 ™ = e* 71 "/ 2 ™, so that = {p}. We have 


7 ™,!— y: 


j=o (2 m l e j — Pi) 2 


- |AE 2m ,2(p)| + 


(2 ,n |p — p|) 2 


= 1 + 


|2 m sin( 7 r/ 2 m )| 2 


where we have used (I2.12p and used that p is a 2 m -root of —1. Since 2 m sin( 7 r/ 2 m ) > 2 then 
7 m. 1 E 5/4. For the other values of / = 2,..., 2 m_1 , we have 


• If / = 2 N, we have for any j > 0 that w Fml (e 2j+ i) = w Fml (e 2j ) = WE 2rn _ h2rn _ 1+N (ej). 
Pairing the indices in (12.24)1 as 2 j and 2 j + 1 with j = 0,..., 2 n_1 + N — 1, we deduce 


7m, l 


7m— 1,JV 


• If l = 2N + 1 with N > 1, we may write 


7m, l 


1 

4m-1 


K -1 


E 


|ejf — ej\ 
w F m , l+ i(ej )| 2 


< 



1 / ^ 


j=0 


jeAT — Cj\ 

w F m , i+ 1 (e 7 ) | 2 


7m—l,JV+lj 


where we have again paired the indices by 2j and 2 j + 1 for j = 0,..., 2 n + (N + 1) — 1 
and used e 2 j+i = — e 2 j and the identity \a + 6| 2 + \a — b\ 2 = 4 for any a, b G dU. 
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Therefore 


7 m,l < 1 <m, 1 < l < 2 m \ 

where (a m ;) i<m is the sequence that saturates the previous inequalities and hence is 
dehned by the following recursion: 

n > 1, N — 1,..., 2 m_2 , 
n> 1,N = 1,..., 2 m ~ 2 — 1. 

The sequence (a m> z) has no dependance on m and it is equal, in the sense a m> i = ai, to the 
sequence (a;);>i which satisfies the recursion: a\ = 1, a, 2 N = ajv/2, a 2 Ar+i = ajv+i- Since 
0i(l) +p(l) = 1, <Ji(2N) + p(2N) = ui(N) +p(N) + 1 and 

0i (2 N + 1) + p(2iV + 1) = 0i (2iV + 1) = 0i (N) + 1 = 0i (N + 1) + p(N + 1), 

then an immediate induction shows that ai = which finishes the proof. ■ 

3 3^-Leja sequences on [—1,1] 

3ft-Leja sequences were introduced and studied in [3], Such sequences are simply dehned 
as the projection, element-wise but without repetition, into [-1,1] of Leja sequences on U 
initiated at 1. More precisely, given E = (ej)j >o a Leja sequence on U initiated at 1, the 
!R-Leja sequence R = 0 associated with E is obtained progressively by: r 0 = 9ft(eo) = 1, 

J(0) = 0 and 

r k = 3where J(k) = min{j > J{k — 1) : 3?(ej) ^ Rk}-, k > 1. (3.1) 

This means one projects e 3 if and only if e 3 ^ e7 for all i < j. The projection rule that 

prevents the repetition is provided in [4j Theorem 2.4], One has 

OO 

R — 5R(5), with S := (1, —1) A A Eyjpi+V- 1 - (3-2) 

3 = 1 

Using a simple cardinality argument, see [4J Theorem 2.4] or [7J Formula 40], this implies 
that the function J used in (13.ip is given by: J(0) = 0, J(l) = 1 and 

J(k) = 2 n + k - 1, n > 0, 2 n + 1 < k < 2 n+l + 1. (3.3) 

In view of (13.21) and the properties of Leja sequences on U , any JULeja sequence R satisfies 
r 0 = 1, ri = —1, r 2 = 0 and r 2 j- \ = — r 2 j for any j > 2. An accessible example of an 3?- 
Leja sequence is the one associated with the simple Leja sequence given by the bit-reversal 


a m ,i = 1, m > 1 and 


® m,2N ®m- l,Jv/2 

0"m,2N+l = 


13 



enumeration (12.3p . We have shown in [6] that R = (cos(0j)) j>o where the sequence of angles 
(0fc)fc>o is defined recursively by 0o = 0, 0i = 7r, (p 2 = 7t/ 2 and 


4>2j~l — 02 j — 4>2j-l + T 1 " 5 j > 2. (3.4) 

This recursion provides a simple process to compute an 3^-Leja sequence. We can also 
construct a Leja sequence by simply using the recursion r 0 = 1, rq = — 1, r 2 = 0 and 


r 2j-i = 



r 2 j = -r 2 j- i, J > 2. 


(3.5) 


One can check that the last sequence is obtained from the Leja sequence F which is con¬ 
structed recursively by F\ = {1} and F 2 n +1 = F 2 n Ae^F 2 n. Both 3?-Leja sequences R satisfies 
2rg — 1 = 1, 2 r 2 — 1 = —1 and more generally 2r^- — 1 = 2r|-_ 1 — 1 = rj for any j > 2, 
thanks to the trigonometric identity 2 cos 2 (0/2) — 1 = cos(0). This shows that in both cases 
R satisfies the property 


R 2 = R where R 2 := (2r%j — 1 )j>o, (3.6) 

In general, given a Leja sequence E in U initiated at 1 and R the associated 3?-Leja 
sequence, we have that R 2 is an 3ft-Leja sequence and it is associated with E 2 which, in view 
of Theorem 12.21 is also a Leja sequence initiated at 1. This result is given in [7J Lemma 3.4] 
and it has many useful implications that we have exploited in order to prove that D k{R) 
grows at worse quadratically. 

For all Leja sequences E on U initiated at 1, the section E 2 n+ 1 is equal in the set sense 
with the set of 2 n+1 -roots of unity, therefore for all 9ft-Leja sequences R , the section /?, 2 « +1 is 
equal to the set of Gauss-Lobatto abscissas of order 2 n , i.e. 

R 2n+1 = { cos(j7r/2”) : j = 0,..., 2 n |, (3.7) 

in the set sense. This set of abscissas is optimal as far as Lebesgue constant is concerned, in 
the sense L^ 2 „ +1 ~ 21 °g(2 ±10 More precisely, we have the bound 

W +1 <l + -log(2 re ), (3.8) 

7T 

see [13j Formulas 5 and 13]. This suggests that the sequence R might have a moderate 
growth of the Lebesgue constant of its section R k . 

In the paper [4], it has been proved that L fifc = 0(k 3 log(fc)). We have improved this 
bound in EE] and showed that < 8\/2 k 2 for any k > 2. Here we again exploit 
our approach of [0] which, using simple calculatory arguments, relate the analysis of the 
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Lebesgue function associated with Rk to that of the Lebesgue function associated with the 
smaller Leja section that yields Rk by projection. This approach allows us to circumvent 
cumbersome real trigonometric functions which arise in the study A# fc , see EflE], and to take 
full benefit from the machinery developed for Leja sequence on li. 

Remark 3.1 Without loss of generality, we assume for the remainder of this section that E 
is the simple Leja sequence in (12.'ll) and R the associated IR-Leja sequence. All our arguments 
hold in the more general case, the assumption is essentially for notational clearness. R allows 
us, in view of (12.51) . to use E instead for E 2 and more generally instead of E 2P which is defined 
by E 2P := ((e 2P j) 2P ) j > o- 

The bound (I3.8[) is sharp and we are only interested in bounding h Ek when A: — 1 is not 
a power of 2. For the remainder of this section, we use the notation 

n > 0, 2 n < k- 1 < 2 n+1 , 0 < l := k- (2 n + 1) < 2” 

K := 2 n+1 + 1, G k :=E K , F K :=E 2n+ p K . 

We should note that in [7] we have used k' and F k to denote l and F K . In view of ( 13.31) . 
we have K = J(k), so that E K is the smallest section that yields Rk by projection into 
[—1,1]. We denote by L 0 , L i, L 2 , - - • , Lk -i the Lagrange polynomials associated with Ek- 
The inspection of the the proof of |7J Lemma 6] shows that for z G dU and x = 9?(z), 

K—l I / \ I 

Ai4<r) < Ak(z) +Ak(z), 7 k(z) ■■= \w Fk (z)\ J j _ ■ (3.10) 

j = o \ W Fk \ e j)\ 

In the proof of [7J Lemma 6], we have bounded the functions |t/;^J/|iup fc (ej)| in the previous 
sum by 2 n+ h~R l ), This implied the result of [7, Theorem 5], namely L^ fc < 2 n+ i~R l)l F EK . 
In view of the new bound (12.231) and the facts that p(K) = p(l), o\ (K) = 1 + <j\(l) and 
K = 2" + fc - 1 < 3 x 2", the previous bound implies 

,— 3n—3p(Z)+tr 1 (Z) 

L. Rk < 12v3 2-3- L 2P (o, k> 1, (3.11) 

where L 2 p is bounded as in (I2.22j) . We propose to improve slightly the previous inequality 
by applying rather Cauchy Schwartz inequality when bounding the function 

Theorem 3.2 Let R be an IR-Leja sequence and n, k and l as in (13.91) . We have 

L/4 < 6\/5 2"+‘*( j )-p(1)l 2P ( 0j (3.12) 

where L 2 p is bounded as in (12.22)) . 
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Proof: In order to lighten the notation, we use the shorthand p in order to denote p(l). We 
introduce l' and K' and Fk> dehned by 


V := 1/ 2 P , K' ■= K/2 P = 2 n ~ p+1 + V, F K < := E 2 n- P+ i >K ,. 


The sequence E satishes E 2 = E and one can check that wp K (z) = wf k ,{z 2P ). Also by 
e 2 j = e lj+i = e ji one has ( e 2 pj+q) 2P — e j for any q = 0 ,..., 2 P — 1 . Moreover, if Mi,..., Mk >-i 
are the Lagrange polynomials associated with E K <. then 

2 P —1 

\L 2 vj+q(z)\ < L 2 P Mj(z 2 ), j — 0,... , K' — 1, 

g=0 

see the proof of [ 6 ], Theorem 2.8]. Therefore by pairing the indices in the sum giving 7 k by 
2 p j + q for j — 0,..., K' — 1 and q = 0,..., 2 P — 1, we infer 


K.' 

1k{z) < (W (01 ^2 I L ) L 2P = L 2P 7 K’iOl with ^ = - 2P - 

In view of (I3.10L this implies that L R k < 2L 2 p sup c &u 7ir'(0• Applying Cauchy Schwatrz 
inequality to 7 k> and using that Fk> is an Z'-Leja sequence, we have for any £ e U 


TWO < 1 ) 1/2 (X] |M,(0| 2 ) 1/2 = 2"‘ (, '> + ”-'' v ^TT7 A e ,„ 2 (0. 


K'-l 


i =0 \ w F K >( e j)\ 


3=0 


where 7 n - P +i,i' is dehned as in (12.241) with m = n—p + 1 and A e k ,,2 is the quadratic Lebesgue 
function associated with E hr view of the bounds we have for these quantities and in view 
of = 1 + cri(Z') and <J\(V) = Ci(Z), we get 


7 *-(£) < 2 n(,>+n ~ P 2^ 


1 < 3\/5 2 ai{l)+n ~ p . 


The proof is then complete. 


_ _ <t\ (l)-\-p(l)—n 

The bound in (13. 12ft improves the bound in (13. lip by 2 2 . The bound can also 

yield linear estimates for L^, for instance when l is such that 2 °W) - M)L 2 p (z ) < 1 , which is 
the case if for example p(l) > 2a\(l). However, if 0 < l < 2 n is the integer with the most 
number of ones in the binary expansion, i.e. oy(/) = n or l = 2 n — 1 and k = 2 n+1 , we merely 
get the quadratic bound 

L r4 < 6^5 2 2n = ^ k 2 . (3.13) 

In |3], section 3.4, it is shown that for the values k = 2 n , in other words Rk is the set of 
Gauss-Lobatto abscissas (j3.7[) missing one abscissa, one has L R k > A R k (rk) = k — 1. As a 
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consequence, the growth of 1L Rk for k > 1 can not be slower than k. However, for this case, 
we can prove < 3k, see (14.lip , showing that (13.131) is rather pessimistic. 

The estimate in (13.12)1 is logarithmic for many values of the integer k. For instance, if 
k = ( 2 n + 1 ) + 2 n ~ p k' for some p = 1 ,..., n and some k' = 0 ,..., 2 P — 1 , then we have 
l = 2 n ~ v k!, so that n — p < p(l ) < n and <7\(l) = cri(k') < p implying that 

h Rk < 6\/5 2 2p L 2 P{ i) < 6\/5 2 2p L 2 . < 6 a/5 2 2 p ^(log(2 n ) + 9/4). (3.14) 

For a small value of p, the previous estimate is as good as the optimal logarithmic estimate 
2lo ^- k ) for large values of n. Given then p fixed, one has 2 P intermediate values between 2 n +1 
and 2 n+1 + 1 , which are the numbers k = ( 2 n + 1 ) + 2 n ~ p k' for k' = 0 ,..., 2 P — 1 , for which 
the Lebesgue constant is logarithmic. This observation can be used in order to modify the 
doubling rule with Clemshaw-Curtis abscissas in the framework of sparse grids, see HU- 

4 Growth of the norms of the difference operators 

In this section, we discuss the growth of the norms of the difference operators Ao = Iz t and 
A k = Iz k+1 ~ Iz k for k > 1, associated with interpolation on Leja or 3ft-Leja sequences. We 
are interested in estimating their norms O*, defined in (1 1. 6 j) . Elementary arguments, see |7j, 
show that 

D k (Z) = (l + A Zk (z k )) sup j WZk ^l , k> 1 . (4.1) 

v / zex \wz k [Zk)\ 

In particular if Z is a Leja sequence on the compact A", then 

Dfc(Z) = 1 + A z k (zk)- (4-2) 

In [ 6 ], we have established that \E k {&k) < k if E is a Leja sequence on U initiated at dU , 
which implies D k{E) < 1 + k. Here, we improve slightly this bound. As for the improvement 
of (12.171) into (12.23]) . we have 

Theorem 4.1 Let E be a Leja sequence on the unit disk initiated at eo € dU , One has 
Oo(E) = 1 and 

Vk(E)< l + /A y ( 2 „W-i) L 2r(w ( 4 , 3 ) 

For 3fJ-Leja sequences R on [—1,1], we have shown in [7] using a recursion argument based 
on the fact that R 2 defined as in (13.61) is also an 3ft-Leja sequence, that 

D fe (i?) < (1 + k) 2 , k> 0. (4.4) 
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In view of the new bounds obtained in this paper for Lebesgue constant of 9?-Leja sections, 
the previous bound is not sharp. Indeed, we have D*, < + L^_i < 12\/5 k 3 ^ 2 , for k 

such that l = k — (2 n + 1) < 2 n / 2 . We give here a sharper bound for D k{R)- We recall 
that up to a rearrangement in the formula (14.11) . see [7J for justification, we may write 
the quantities O k(R) in a more convenient form for JGLeja sequences. We introduce the 
polynomial W Rk := 2 k w Rk , we have 

D t (fi) = 2fo(R) sup \W Rt (x)\, 0 t (R) := 4±^sM , ( 4 , 5 ) 

We have already proved in [71 Lemma 7] that 

p 2 n (R) = 1/4 and /3 fc (i?) < 2 ffo(fc) - p(fc) - 1 , for k ± 2 n . (4.6) 

Here we provide a sharper bound for D k(R) by slightly improving the estimate zjyiW+HW 1 
that we have established in m for sup iehl l] \W Rk (x)\. 

Lemma 4.2 Let R be an tR-Leja sequence in [—1,1], n > 1, 2 n + 1 < k < 2 n+1 + 1 and 
l = k — (2 n + 1 ). One has sup xe [_ 1;1 ] \W Rk (x)\ < 2 n+3 if k = 2 n+1 , else 

sup \W Rk {x)\<2 2ai{ - k)+p(k) - 1 . (4.7) 

®e[—i,i] 

Proof: We use the notation K, Cf,. and F K as in (I3.9[) and introduce G'/,. + 1 := E K+l and 
Fk +i := E 2 n+i > K+i- In view of [7, Lemma 5], one has for z G dU and x = SR(z) 

\W Rk (x)\ = \z 2 - l\\w Gk {z)\\w FK (z)\ = \z -z\\w Gk (z)\\w FK (z)\. 

Also since \z — z\ < \z — ex\ + \z — e^|, then 

\W Rk (x)\ < \w Gk+1 (z)\\w FK (z)\ + \w Gk {z)\\w FK+1 {z)\. 

In the two previous inequalities, one has F K = 0 and uj Fk = 1 in the case k = 2 n + 1. We 
have that Gk, Gk+ 1 , Fk and Fk+i are all Leja sections with length K, K + 1, l and l + 1 
respectively. Therefore, by the second property in Theorem 12.21 

\W Rk (x)\ < m i n ^2 1 +CTl(K ) +f 7 i ( 0 ,2 <Tl( ^ + 1 ) +CT i(') + 2 <Tl W+<n 6 +i)^ = 2 2 + f 7 i(b m i n (2 <T dO ) 2 CTl ^+ 1 )), 

where we have used <j\(K) = 1 + &i(l) and o\ (K + 1) = 1 + < 7 i(/ + 1) since K = 2 n+1 + l 
and l < 2 n . If k — 2 n+1 , i.e. I = 2 n — 1, then \W Rk (x)\ < 2 3+n . Else by k = 2 n + (Z + 1) and 
0 < l < 2 n - 1 , 

tj\(k) — 1 — (?\[l T 1) and cri(k) — 2 + p(k) — a\{k — 1) — 1 = 
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Therefore 

\W Rk {x)\ < 2 2fTl(fc)+p(fc) - 1 min(2' 1+p(fc) ,l), 

which completes the proof. ■ 

By injecting the estimate of the previous lemma and the estimate of (j4.6[) in formula 
(USD and by using the identity a 0 (k) + <i\(k) = n + 1 for 2 n < k < 2" +1 , we are able to 
conclude the following result. 

Corollary 4.3 Let R be an R-Leja sequence in [—1,1]. The norms of the difference operators 
satisfy, D 0 (R) = 1 and for 2 n < k < 2 n+1 

B k (R) < 2 ai{ - k) 2 n (4.8) 

The previous estimates can be used in order to provide estimates for L Rk that can be 
sharper than (13.121) . We have A k = I k ~ R-i, therefore 

|Lhj. +1 — LrJ < D k{R), k > 1. (4.9) 

In particular, the estimate in the previous corollary combined with the sharp bound (13.81) 
implies that for the value k — 2 n , we get 

h R . < 1 + - log(2 n ) + 2 n+1 < 3 k (4.10) 

7r 

This shows that in the case k = 2 n which corresponds to R k being the set of Gauss-Lobatto 
abscissas (|3.7[) missing one abscissas and for which h Rk > k, the previous bound is satisfac¬ 
tory. This also confirm that the estimates (13.121) is indeed pessimistic in this case, see the 
inequality (13.131) . This added to the observed growth of for values k < 128, Figure l4~Tl 
suggests that the bound 

Ui k <3 k, k > 1, (4.11) 

might be valid for any 9?-Leja sequence R. We conjecture its validity. 

In Figure 14.11 we also represent for the values k < 128, the growth of the Lebesgue 
constant L Rk (in blue) and the estimate \/k{2 a d k ) — 1) (in red) which multiplied by 3 bounds 
L E k , see (12 .1 7j) . We observe that the regular patterns in the graph of k which reveals 

the particular role of divisibility by powers of 2 in k, is caught by the estimate. The worst 
values of appear for the values k = 2 n — 1 for which it was proved in [13] that L E k = k 
and which is also equal to y/k[2 <T d fc ) — 1) since <J\(k) = n. 
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Figure 4.1: Exact Lebesgue constants associated to the ^-sections of the Leja sequence E 
and the assciated 9ft-Leja sequence R for k — 1, 3,..., 129. 
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